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ABSTRACT 

We report on the present stage of SN 1987A as observed by the Chandra X-ray Observatory. We 
reanalyze published Chandra observations and add three more epochs of Chandra data to get a consis- 
tent picture of the evolution of the X-ray fluxes in several energy bands. We discuss the implications 
of several calibration issues for Chandra data. Using the most recent Chandra calibration files, we 
find that the 0.5-2.0 keV band fluxes of SN 1987A have increased by ~ 6 xl0~ 13 erg s _1 cm~ 2 per 
year since 2009. This is in contrast with our previous result that the 0.5-2.0 keV light curve showed a 
sudden flattening in 2009. Based on our new analysis, we conclude that the forward shock is still in 
full interaction with the equatorial ring. 

Subject headings: ISM: supernova remnants - supernovae: individual (SN 1987A) - X-rays: individual 
(SN 1987A) - radiation mechanisms: thermal 



1. INTRODUCTION 

Stellar winds of massive stars leave their imprints on 
the surrounding medium before they explode. After the 
explosion, the outer shocks of the supernova probe the 
circumstellar medium, thereby mapping the late stages 
of the stars evolution. 

Because of its youth, the remnant of SN 1987A pro- 
vides a unique opportunity to investigate these late 
stages. The case of SN 1987A is especially interesting, 
as it exhibits the enigmatic three ring structure visible 



in optical images (Crotts et al. 19891. This structure 



probably originates from the interaction of a wind from 
the slow red supergiant p hase with the faster wind from 
the blue super giant phase ( Luo &: McCray|1991 Wang &; 
Mazzali 1992). The morphology ot this structure might 
pro vide clues to a binary-merger scenario for the progeni 



tor (Blond jnfc Lundqvist||1993 Morris & Podsiadlowski 
20071 I2OO9I). However, this three ring structure might 



20U9[) . However, this three ring structure might 
also result f rom mass-loss from a fast-rotating star ( Chita 



et al.|[2008 1 . Currently, the explosion is sweeping up the 
inner equatorial ring that was formed by the late stages 
of the star's evolution. 

As the inner equatorial ring is ab out a hundred time s 
denser than its surroundings (e.g., Mattila et al.| 2010), 
the interaction of the explosion with this structure has a 
direct impact on the evolution and observational proper- 
ties of the remnant. In 1995, the outer ejecta began to 
light up this structure. By October 2000, 12 spots had 
lit up at visible wavelengths ( Sugerman et^|2002| . In- 
frared observations indicate that t he lull in teraction with 
the ring started around mid-2002 ( Dwek et al.|2010[ ). A 
year and a half lat er, the expansion in X-rays slowed 
down dramatically (Racusin et al. 2009) while, at the 



same time ; the soft X-ray light curve turned up ( Park 
eTaLl[2005 ). 



Base d on past Chandra observations, Zh ekov et al. 
(2005) established that the X-rays are emitted from a 
flat spatial structure. A model consisting of three com- 
ponents has been proposed for interpreting the high reso- 
lution X-ray spectroscopy from SN 1987A (Zhekov et al. 
2010 Dewey et al. 2012). The first component is emit- 



ted by material in the dense ring that has been heated 
by a slow shock (500 km s~ x ). This component has a 
low electron temperature (~ 0.3 — 0.5 keV). The sec- 
ond component is emitted by the shock-heated circum- 
stellar medium located above and below the equatorial 
ring. These shocks have not been slowed down as much 
and therefore the plasma is hotter. Also, the line emis- 
sion of this component has to be broad (t housands of 
km s~ x ). This allowed Dewey et al. (20121 to estimate 
the current contribution of this component to be of the 
order of ~ 20% of the total emission in the 0.5-2.0 keV 
band, based on high resolution High Energy Transmis- 
sion Grating (HETG) spectra and Reflection Grating 
Spectrometer (RGS) spectra. Further investigation of 
the X-ray spectra of SN 1987A shows that there is need 
for a third component with a higher (~ 3 keV) temper- 
ature, but without broad line emission. Zhekov et al. 
(2006) proposed that this is shocked plasma, reheated 



by a reflected shock caused b y the interacti o n wit h the 
equatorial ring. Alternatively, Dewey et al. (20121 pro- 



posed that this high temperature component arises from 
a fast shock, traveling through the uniform equatorial 
ring material (in contrast to the low temperature com- 
ponent, which they consider being emitted by shocked 
denser clumps in the equatorial ring). 

The Chandra X-ray Observatory has observed SN 
1987A at a regular pace of about two observations per 
year, resulting in a denser coverage of the evolution of the 
remnant th an with any other X-ray observatory. Park 



et al. (2011) summarized the observations through 2010. 
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In the 0.5-2.0 keV light curve, they found a break (i.e., 
the pace with which the flux increases suddenly slowed 
down) around 8200 days after the explosion (~ Y2009). 
Based on this result, they concluded that the shocks re- 
sponsible for the soft X-ray emission must have recently 
encountered less dense material. 

However, the May 2012 revision to the ACIS calibra- 
tion (CALDB 4.4.100 introduced a significant change to 
the time-dependent model for the transmission of the op- 
tical blocking filter. The new model predicts lower trans- 
mission than the previous model (a 15% change in the 
optical depth in 2012), reflecting an apparent increase in 
the rate of deposition on the filter. 

In this paper, we rec onsider the light curve break re- 
sult (Park et al. 2011) by reanalyzing previously pub- 
lished epochs using the newest calibration and by adding 
three recent observations (2011 March and September, 
2012 A pril) . We show that the break seen by Park et al. 
(2011 1 is due primarily to the (unmodeled) increased rate 
of deposition on the optical blocking filter and not to a 
different phase in the evolution of the remnant. 

2. DATA AND RESULTS 

Table JT] lists the 23 monitoring observations used for 
this study We chose not to use the observations earlier in 
the mission (Obslds 1387 and 122 at 4608 and 4711 days 
after the explosion), as these observations were taken at 
a focal plane temperature of -110° C and no charge trans- 
fer inefficiency adjustments are available for observations 
taken at this temperature^] However, we report the spa- 
tial distr ibution of the X-ray data at these epochs in sec- 
tion 2.1 Since the remnant has thus far brightened by a 
factor of ~ 34 during the Chandra mission, the observ- 
ing configuration was changed several times to mitigate 
photon pileup effect^] First, this was done by read- 
ing out an increasingly smaller subarray of the ACIS- 
S3 chip. All observations from day 7799 onward used 
the ACIS-HETG configuratiorj^J which produces both a 
zeroth-order image and a dispersed spectrum. The dis- 
persed spectrum provides an invaluable cross-check for 
our ACIS results. From day 8434, the observation was 
moved to the bottom of the chip, to decrease the frame 
time to Is. The time of this transition coincides with the 
accelerated pace of the increase of the contamination on 
the filter. For the data reduction, we largely follow the 
met ho d s that were used before on Chandra data of SN 
1987A (IBurrows et al.||2000[ IRacusin et al.||2009[ [Park 



much spatial information from the observations as pos- 
sible, we adopt the same approach as in previous studies 



( Burrows et al.||2000| |Racusin et al ||2009[ |Park et al 
2Ull|. We remove the pixel randomization that is usu- 



ally added by the Chandra software. Additionally, we 
improved the spati al resolution using split-pixel events 



(Mori et al. 2001). We then deconvolve the image 



for the point spread function of the telescop^] using 
the Lucy-R i chardson iterativ e deconvolution algorithm 
( [Lilcyl[l974j |Richardson|1972[ ). 

For the September 2U11 and March 2012 observations, 
we analyze only the longest exposure, to avoid distortions 
caused by residual alignment errors between exposures. 
Figure [T] shows the deconvolved images of all epochs of 
Chandra observations to date. To characterize these im- 
ages, we fit the i mage with a model that consists of four 
lobes and a ring (Racusin et al. 2009 1. Figure [2] shows the 
best-fit radii of all the images considered in this study. 

2.2. Spectra 

Spectra were extracted from each observation within a 
circular aperture (radius 4'.'38) sized to encompass nearly 
all the light from the target in al l epochs, using CIAO 
version 4.4 (Fruscione et al. 2006) and CALDB version 
4.4.10. Background spectra were extracted from each ob- 
servation within an annular region (radius 6'.'2 to 12'.'4). 

Even though we adjusted the observing modes to min- 
imize pileup, PIMMS suggests that a modest pileup is 
present in many of the observations. We first tried to 
take this into account during the fi tting process by us- 



ing the pileup model described by Davis 



(2001). This 



et al.||20TT| 

nored the flags set by the Chandra afterglow filter (sta- 
tus=000000000000xxxx0000000000000000) as this filter 
tends to filter r eal events for bright sources (e.g. 
ley et al.|[2003| ). 



approach adds to the overall model several pileup pa- 
rameters which then have to be either assigned appro- 
priate values or included in the fitting: for details see 
the "Chandra ABC Guide to Pileup'fl As an exam- 
ple, the parameter nregions encodes the spatial size of 
the source in units of 3 x 3 pixel cells. Based on the 
ACIS image we can estimate that nregions should be in 
the range of 2 to 5, but we cannot determine the exact 
value a priori. Including nregions in the fitting how- 
ever adds a degeneracy between nregions and the other 
model parameters such that the fit and its flux are not 
well constrained. 

An alternative, parameter-free strategy for addressing 
pil eup in the analysis of ACIS data has been described 
by Broos et al. ( |2011| Appendix A). First, a high-fidelity 



We filtered tor background flares, and ig- simulator of the ACIS detector (Townsley et al. 2002a 



includi ng charge transfer inefficiency (Townsley et alj 



Towns- 



2.1. Imaging 

With an angular diameter of 1'.'5, SN 1987A is barely 
resolved in raw Chandra ACIS images. To extract as 

1 http : //exc . harvard . edu/caldb/downloads/Release_notes/ 
\CALDB _v4 .4.15. htmT|~ 

z http : //exc .harvard. edu/ ciao4. 4/why/cti .html 



6 http : //exc .harvard. edu/ciao/ahelp/acis_pileup. html 

4 Use of the HE'l'G is an excellent pileup mitigation strategy, 
since the count rate in the zeroth order is reduced by a factor of 9. 



[2002bJ, is used to infer an input photon spectrum that 
would cause ACIS to produce the observed spectrum and 
event rate under the specific configuration of the obser- 
vation. That photon spectrum is very flexible; no spec- 
tral model is assumed. The spatial distribution of the 
simulated photons incident on ACIS was produced by 
the MARX packag<{H using the corresponding SN 1987A 
morphology model aescribed in |2.1| MARX and the 
CCD simulation were configured to match the parame- 
ters of the observation, such as the target location on the 

5 Using MARX version 4.4.0 (space.mit.edu/cxc/marx) 

6 http:/ /exc. harvard.edu/ciao/download/doc/pileup_abc.pdf, 
Section 4.1: Correcting Imaging Observations. 

7 http://space.mit.edu/cxc/marx/ 
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TABLE 1 
ACIS-S3 Observations 

Obsld(s) SN age Instrument frame time 0.5-2.0 keV flux" 3.0-8.0 keV flux" 

days [seconds] [10 -13 erg s _1 cm -2 ] [10 -13 erg s _1 cm 



1967 


5036 


ACIS-S3 


3.2 


1044 


5175 


ACIS-S3 


3.2 


2831 


5406 


ACIS-S3 


3.1 


2832 


5560 


ACIS-S3 


3.1 


3829 


5790 


ACIS-S3 


3.1 


3830 


5979 


ACIS-S3 


3.1 


4614 


6157 


ACIS-S3 


3.1 


4615 


6359 


ACIS-S3 


1.5 


5579 & 6178 


6530 


ACIS-S3 


0.4 


5580 & 6345 


6713 


ACIS-S3 


0.4 


6668 


6914 


ACIS-S3 


0.4 


6669 


7094 


ACIS-S3 


0.4 


7636 


7270 


ACIS-S3 


0.4 


7637 


7445 


ACIS-S3 


0.4 


9142 & 9806 


7625 


ACIS-S3 


0.2 


9144 


7799 


HETG 


1.1 


10221 & 10852-10855 


8000 


HETG 


1.1 


10222 & 10926 


8233 


HETG 


1.1 


12125 & 12126 & 11090 


8434 


HETG 


1.0 


13131 & 11091 


8618 


HETG 


1.0 


12539 


8796 


HETG 


1.0 


14344 & 12540 
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HETG 
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9175 


HETG 


1.0 
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(0.85) 
(1.01) 
(0.86) 
(0.71) 
(0.76) 
(0.82) 
(0.83) 
(0.88) 
(0.83) 
(0.94) 
(0.96) 
(0.85) 
(0.89) 
(0.77) 
(0.63) 
(0.92) 
(0.81) 
(0.89) 
(0.90) 
(0.85) 
(0.98) 
(0.88) 
(0.97) 



Observation log of the Chandra S3 observations of SN 1987A. 
corrected fluxes and the non-pileup corrected fluxes. 

focal plane, frame time, and dead area/time[^] 

Second, with the input photon spectrum fixed, the sim- 
ulation is repeated with only a single photon arriving in 
each frame, eliminating pileup. The resulting simulated 
spectrum, which we refer to as "pileup corrected" , is then 
fit using standard methods and calibration files. Figure 

t shows the original and pileup-corrected spectrum of 
e observation that suffered worst from pileup (Obsld 
4614). We applied this method to all ACIS observations 
reported here. The flux correction attributed to pileup in 
each observation is reported in Table[l] This Table shows 
that the correction for the soft band is more smooth than 
the correction of the harder band. This probably indi- 
cates that our pileup correction method is more stable 
for the soft band than the hard band. 

2.2.1. Light curve 

We fitted the resulting pileup-corrected spectra with 
an absorbed two component spectral model, using the 
X-ray spectral fitting package XSPEC, version 12.6.0 
(Arnaud 1996). The spectral model included a compo- 
nent m collisional equilibrium ionization (vequil) and a 
non-equilibrium ionization component (vpshock) using 
neivers 2.0, updated with a tomic data for inner shell 
processes (c.f. |Badenes et al.p006[ ).We fixed the abun- 
dances to those tound m the high resolution X-ray LETG 
spectra investigated by Zhekov et al. ( 2006 ) and we fixe d 
the absorption to 2.35xlO zi cm" 2 (Park et al.||2006|). 



1 Numbers between parentheses are the ratios of the pileup 

Figure 14] shows an example of one of the fits. We mea- 
sured the flux in the 0.5-2.0 and 3.0-8.0 keV energy bands 
from the best-fit model for all spectra and estimated flux 
uncertainties for a 68% confidence level. Figure [3] shows 
the evolution of the fluxes in the low and high energy 
bands for the observations corrected for pileup effects. 
For this study, we do not need the high signal-to-noise 
that would be required for a detailed spectral analysis. 
Therefore, we decided to only correct the longest ob- 
servation per epoch for pileup even if there were more 
observations available in the same epoch. This explains 
the large error bars of the flux at day 7997, as this epoch 
was divided into 5 observations. 
For comparison, we overplotted the fluxe s of the last 



six epochs of XMM Newton EPIC pn fluxes (Maggi et al. 
:). Note that this instrument has been shown to be 



2012) 



extre mely stable (e.g., Sartore et al. 2012 Plucinsky et al. 
20121. 



8 See discussion of F RACEXPO at |http : //exc ■ harvard . edu/| 
ciao4 . 4/why/acisdeadarea. html 



2.2.2. Chandra grating observations 

After the previously described procedures, the 0.5-2.0 
keV light curve increases more or less linearly from day 
6400 on, but there still appears to be a break around day 
8200, which do es not appear in th e XMM-Newton EPIC 
pn light curve ( Maggi et al.||2012 ). 

Since the pileup corrections ( Table fll) after the break 
are relatively low (7-11%) compared to some earlier ob- 
servations, insufficient pileup correction is an unlikely ex- 
planation for the break. 

We can confirm this directly using the data taken with 
the ACIS- HETG configuration, since the dispersed spec- 
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Fig. 1. — Estimated intrinsic morphology (corrected for the Chandra PSF) of SN 1987A at all epochs utilized in this study, proceeding 
from upper left to lower right. The image scaling is linear and the brightness scales with the 0.5-8.0 flux in the observation. The scaling 
limits of the two top rows are offset by a factor of 3 for visibility. In the images, north is up, east is left. 
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Fig. 2. — Radii of a toroidal model (j ]2.1[ l fit to each image in 
Figure 1. 

tra in each of these observations have no pileup, they 
provide invaluable cross-checks for our pileup-corrected 
ACIS fluxes. The dispersed spectra were extracted us- 
ing the standard CIAO tools and extraction sizes. The 
Medium Energy Grating (MEG) plus and minus data 
are combined and the fluxes in the 0.79 to 2.1 keV range 
are determin ed directly from the fl ux-corrected data, as 
was done in |Dewey et al. (2012 1. In addition to the 
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Fig. 3.— Chandra ACIS light curves for the soft (0.5-2.0 keV, 
higher curve) and hard X-ray (3.0-8.0 keV, lower curve) bands, 
with pileup correction (§2.2). XMM-Newton EPIC pn fluxes are 
shown (brown) for comparison. Uncertainties shown are lc. The 
XMM-Newton EPIC pn fluxes are displayed with 3 cr-errors. 



the deep observations taken in 2007 (day 7335) and 2011 
(day 8774). These observations have been described in 
Dewey et al. (2012). The light curve obtained from the 



monitoring observations listed in Table [T] we analyzed 



dispersed spectra (Figure [5]) shows a break around day 
8200 similar to that seen m Figure [3] Since this break 
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Fig. 4. — Original (black) and pileup-corrccted spectrum of Ob- 
sld 4614 with best-fit model overplotted. Note that the piledup 
spectrum is substantially harder than the pileup-corrected spec- 
trum. 
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Fig. 5.— Chandra HETG soft band (0.79-2.1 keV) light curves 
in the dispersed spectra. Note that around day 8200 there is a 
suggestion of a break. The lines are fitted to the data points for 
similar instrument setups individually. Uncertainties shown are 
la. 

occurs at the same time that we moved the aimpoint 
on the CCD, we cannot disentangle any effects caused 
by that change from those caused by the increased con- 
tamination deposition rate, but we can be confident that 
this break is not caused by the increased count rate and 
insufficient pileup correction. 



3.1. 



3. DISCUSSION 
Soft X-ray emission 



Haberl et al. (2006) predicted a flattening of the 0.5- 
2.0 keV light curve around day 7000 and predicted the 
light curve to reach its m aximum around day 8000. Also, 



Borkowski et al. (1997) predicted the soft X-ray light 
curve to reach a maximum around day 9000, but with 
a flux of at up to a factor of ten abo ve the current X- 
ray flux. Indeed, Park et al. (20111 showed that the 



flux in the 0.5-2.0 keV energy band flattened around day 
8200 and argued that this flattening is a result of the 
transmitted shocks entering a lower-density region. 

The improved calibration of the optical blocking filter 
now available removes a substantial part of this flatten- 
ing. However, there still appears to be a break in the 
Chandra light curves, both in the piledup data (bare 
ACIS and zeroth order HETG) and in the pileup-frec 
dispersed data. Note that all the Chandra data pre- 
sented here rely on the same calibration model for the 



optical blocking filter. The XMM-Newton EPIC-pn data 
show no evidence for for a sudden break around day 8200. 
Given that this break is visible in both dispersed and ze- 
roth order, but not in the XMM-Newton EPIC pn data, 
we conclude that this is not likely an effect of an insuffi- 
cient pileup correction, but rather is caused by a miscal- 
ibration of the increase in contamination on the optical 
blocking filters, combined with a change in instrument 
setup. Note that this contamination has a wavelength 
dependent effect on the spectrum and hence has influ- 
ence on the best-fit spectral parameters. Also, this con- 
tamination has less of an effect on the 3.0-8.0 keV flux. 

3.2. Imaging 

Figure [T] shows the evolution of SN 1987A's morphol- 
ogy. To interpret these images, one has to keep in mind 
that some of the structure in the images might be an 
artifact of the deconvolution method we used to correct 
for the extension of the point spread function of Chan- 
dra. Moreover, the point spread function of Chandra is 
not entirely symmetric, as since December 2000, there 
appeared a hook-shaped extension to the point spread 
function^] which is not accounted for in MARX 4.4. This 
hook contains about 6% of the total flux and will proba- 
bly show up after the deconvolution as a small extending 
structure to the remnant. 

The individual spots on the west side of the remnant 
appear to have merged in the last three epochs. One 
question that comes to mind when looking at Figure [2] 
is "how doe s this relate to t he in ner- ring structure as 



observed by |Jakobsen et ahl ( |1991[ )?" . |Ng et al.| ( 12009 
compared radii of X-ray and radio images of SN 1987. 
and found comparable sizes, with radii larger t han the 
radius of the equatorial ring in 2008. However, |Dewey| 



et al. (2012) showed a strong break in the evolution of 
the soft X-ray light curve after the shocks have overcome 
the equatorial ring, and we do not yet see such evolution 
in the 0.5-2.0 keV light curve presented in Figure [3| 

To estimate how far the shocks progressed through the 
equatori al ring, we need t o kno w when the interaction 



started. Sugerman et al. (2002) carried out a study on 
multiple epochs of Hubble Space Telescope (HST) im- 
ages, ranging from 1994 February and 2000 February. 
The first spot lit up by the blast wave is located in the 
northeast quadrant of the remnant and was detected in 
a HST observations taken 2933 days after the explosion. 
Remarkably, this spot is located 0.06" inside o f the main 
ring as visible in optical ( |Sugerman et al.|2 002 ) . The sec- 
ond spot was only found in a HST observation taken 1350 
days later and by day 4999, 12 spots wer e visible in the 
HST observations (Sugerma n et al.| [2002). A lower limit 
to the width of the ring, as visi ble in the HST im age, is of 
the order of 0'.'12 (~ 10 12 km, |Plait et"aT| 1995 1). Fitting 
a broken linear function to the radii of b igure[S]results in 
a velocity of 8500]^^ km s" 1 until day 59OOI230, after 
which the velocity is 1820^ggQ km s _1 . Assuming that 
the shock wave reached the ring at day 4999 when 12 
spots were lit up (hereby ignoring the spots that were 
hit first, as those were protruding inwards), the radius 
has increased by (1.17 ± 0.06) x 10 12 km (-120% of the 
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width reported by Plait et al.|1995 ) since the interaction 
started. We used a Monte Carlo simulation to estimate 
the uncertainties on this value. Note that this distance 
relies on the assumptio n that the 12 spots reported by 
Sugerman et al. l (120021 ) are part of the ring reported by 
Plait et al.| jiyy5B~ 



Dwek et al. 1 2010| ) argued, based on the evolution of 
the infrared emission, that the forward shock only fully 
started to interact with the equatorial ring at day ~5600. 
This appears consiste nt with the increa se in soft X-ray 
flux around d ay 6000 (fPark et al.l|2005[ ). The latter was 
considered by Dewey et al. 12012) to be consistent with 
the lighting up of the optical spots ~ 1000 days earlier, as 
their models show some delay between the shock impact 
of dense material and the onset of X-ray emission. Nev- 
ertheless, if we take day 5600 as conservative assumption 
for the start of the explosion to interact with the ring, 
we estimate the increase of the radius of the X-ray emis- 
sion since then to be 7.36 ±0.05 x 1 11 km (-74% of the 
width reported by [Plait et al.|1995 |. 

It is not straightforward to interpret the increase of the 
radius of the X-ray emission in terms of the 3 components 
as described in the introduction. The velocity is too high 
for the component w it h low (~ 0.3-0 . 5 keV ) temperature 



(Zhekov et al. 2009). Dewey et al. (2012) reported the 



contribution of a high velocity component to be ~ 20 % 
in the HETG spectra, but its contribution to the broad- 
band flux might be higher. Note that they measured 
the spectral lines to be much broader than 1700 km s . 
Therefore, the measured increase of the radius is likely a 
combination of the increase in radius of all three compo- 
nents. The radius of the X-ray emission is identified in 



Dewey et al. (2012) as following the progress of the for- 
ward shock/contact discontinuity as it moves through the 
uniform- wit h- clum ps equatorial ring. Figure 9 (top) of 



Dewey et al. ( 2012 ) shows the forw ard shock/contact dis- 
continuity locations along with the Racusin et al.| ( 2~009| ) 
radii and there is reasonable agreement. The speed of 



speed of 

the forward shock there is of order 1700 km s -1 . 

Current and future HST observations will reveal how 
far the shock has progressed, and how many of the opti- 
cal spots have merged by now. A combined study with 
optical, X-ray and infrared data will reveal a detailed 
picture of the structure of this equatorial ring, leading 



to better constraints on the last stages in the life of the 
progenitor of SN 1987A. 

4. CONCLUSIONS 

We report our results for the three newest epochs of 
Chandra data on SN 1987A. Together with our reanaly- 
sis of older epochs, we paint a picture of the evolution of 
SN 1987A as witnessed by Chandra in the past 10 years. 
Based on this analysis, we reach the following conclu- 
sions: 

- The increase of the radius of the ring visible in X- 
rays since the explosion first interacted with the 
preexisting equatorial ring is 74-120% of the lower 
limit to the thicknes s of this preexistin g equatorial 
ring as measured by Plait et al. ( 1995 1 . 



Given that the shocks have now traveled through 
a substantial amount of the preexisting equatorial 
ring, a multiwavelength study including high spa- 
tial resolution optical images is necessary to inves- 
tigate how far the shock has traveled through the 
ring. 

The sudden break in the 0.5- 2.0 keV light curve 



around day 8200 reported by Park et al. (20111 



was probably an instrumental effect. The absence 
of a break in the light curve likely also indicates 
that the shocks traveling through the dense parts 
of the ring have not fully overcome the extent of 
the ring yet. 
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